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Abstract 

Universal low-energy properties are studied for three identical bosons confined in two dimensions. 
The short-range pair-wise interaction in the low-energy limit is described by means of the boundary 
condition model. The wave function is expanded in a set of eigenfunctions on the hypersphere and 
the system of hyper-radial equations is used to obtain analytical and numerical results. Within 
the framework of this method, exact analytical expressions are derived for the eigenpotentials 
and the coupling terms of hyper-radial equations. The derivation of the coupling terms is generally 
applicable to a variety of three-body problems provided the interaction is described by the boundary 
condition model. The asymptotic form of the total wave function at a small and a large hyper-radius 
p is studied and the universal logarithmic dependence ~ ln 3 /3 in the vicinity of the triple-collision 
point is derived. Precise three-body binding energies and the 2 + 1 scattering length are calculated. 

PACS numbers: 03.65. Ge, 21.45,+v, 34.50.-s 
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I. INTRODUCTION 




. An additional motivation is roused by specific features of quantum 
llOi l . The experiments with ultra-cold gases in the 2D and quasi-2D 

la. 



Dynamics of few particles confined in two dimensions (2D) is of interest in connection 
with numerous investigations ranging from ultra-cold gases [l , Q, 0, Ul to atoms absorbed 
on a surface 
systems in 2D 

traps have been recently realized 

The description of elementary processes in ultra-cold gases has been attracting great 
interest in the last years and many aspects of the low-energy few-body dynamics in three 
dimensions (3D) have been thoroughly investigated. As a particular important example, one 
could mention the studies of the three-body recombination for spinless bosons 11, I2.li3l.li4j]. 
two-component fermions Il5ll. and particles with internal degrees of freedom in the presence 
of a Feshbach resonance jld Il7 |. Low-energy few-body dynamics in low dimensions is 
less investigated despite the extensive studies. Besides, new phenomena and additional 
complications arise in quasi- low-dimensional geometry if the effect of motion in the transverse 
directions is taken into account (see, e. g. 



18, 



19 







2l|). Among different aspects of 



2D systems, one should mention treatment of the three-body energy spectra 



low-energy scattering of an atom off a dimer molecule [251] . and low-energy three-to-three 



scattering 



26 



22 



23 



2i, 



271 ] . The precise binding energy of four bosons was calculated in [28] and the 



universal law for the iV-boson ground-state energy was discussed in [29]. 

Concerning other 2D and quasi-2D problems, considerable efforts have long been devoted 
to investigation of atoms adsorbed on a surface, including helium atoms on graphite .3c| and 



of a quasi-condensate , measurement of the three-body recombination rate |7( and a vast 



hydrogen atoms on a helium film . In this respect, one should mention observation 



n 



number of theoretical papers, e. g., [lCt l31l 

In the low-energy limit, which is of interest both for practical applications and from 
the general point of view, the description of the few-body system becomes universal, i. e., 
essentially independent of the details of the two-body interactions. Among different results 
of this sort, notice the recent analytical derivation of the universal constants for the zero- 

14, 3j|. The description becomes parameterless if only 



energy three-boson scattering in 3D 



the parameter describin g th e two-body interactions, e. g., the two-body scattering length a, 
is chosen ClS db scale In comparison with the universal description in 3D, it is 
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of importance that the solutions of the three-body problem in 2D , contrary to the 3D one, 
remain regular near the triple-collision point even in the zero-range- interaction limit. This 
regularity implies, inter alia, absence of the Thomas and Efimov effect, which was noticed 
m Q<0] Thus, there is no additional regularization parameter, which was introduced in 3D, 
and the three-body properties in 2D are completely determined by the two-body input that 
provides a completely universal and parameterless description in the low-energy limit. In 
particular, the trimer binding energies and the 2 + 1 scattering length become the universal 
constants, which must be determined with a good accuracy. 

The universal limit corresponds to the limit of the vanishing interaction range r$ so that 
r must be much smaller than any length scale in the system, i. e., the binding energies do 
not exceed the characteristic energy — ^. Practically, the universal limit can be realized by 
adjusting the parameters of interaction to diminish the two-body binding energy, e. g., by 
tuning the position of the Feshbach resonance. To approach the universal limit, one could 
use the dependence on the particle mass, e. g., an interesting way is to study the isotopic 
effect for 2D helium atoms (3]. Generally, the universal limit in 2D appears for a very 
weak potential with a non-positive average as it is known that in this case the two-b ody 
binding energy in 2D becomes exponentially small with decreasing potential strength (37 1. 
It is worthwhile to mention a qualitative difference of the universal properties in 2D and 3D, 
namely, one expects that the both three-body bound states arise simultaneously with the 
two-body bound state, i. e., at the infinite two-body scattering length. This conjecture is 



supported by the calculations 22, 38], which considered the dependence of the three-body 
binding energy on the particle mass and potential strength in the limit of the vanishing 
two-body binding energy. Recall that in 3D an infinite number of three-body bound states 
arise with increasing potential strength at finite values of the scattering length before the 
two-body bound state arises. 

In the present paper, universal properties of three identical 2D bosons are studied within 
the framework of the method which makes use of the boundary condition model (BCM) P, 
I39I for the s-wave inter-particle interactions. The wave function is expanded in a set of 
eigenfunctions on the hypersphere and a resulting system of hyper-radial equations (HREs) 
is used to conveniently treat both the boundary and scattering problem. One of the principal 
advantages is that for the two-body interaction given by the BCM the eigenpotentials of 
HREs are solutions of a simple eigenvalue equation. In addition, the aim of the present paper 
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is to derive all the terms of HREs in the analytical form, thus determining the coupling terms 
via the eigenpotentials and their derivatives over the hyper-radius. Both the derivation and 
the expressions found are generally applicable to a variety of problems; as a matter of fact, 
one can obtain the coupling terms in the analytical form for three particles of arbitrary 
permutation symmetry, with arbitrary masses and scattering lengths, and independently 
of the configuration-space dimension. In particular, analytical expressions for the coupling 
terms of the same form as in the present paper were derived for three identical bosons in 
three dimensions 3^|. The exact expressions are used to analyze the asymptotic behaviour of 
the coefficients of HREs, their solutions, and the total wave function both for large and small 
inter-particle separation. As an important example, the universal dependence of the total 
wave function in the vicinity of the triple-collision point is found. In addition, the explicitly 
known dependence on the channel number is helpful to study the role of the channel coupling 
and to estimate convergence of the results with increasing number of HREs. 

Until now, the numerical calculations of the universal constants have included the early 
calculation j8j of the binding energies of three 2D bosons by solving the momentum-space 
integral e quat ions. Much better precision was obtained by solution of the hyper-radial 
equations [2 3! l^ . which results were not in complete agreement with those of 0. Up to 
date, in the universal limit of zero-range interactions the most precise binding energies have 
been found by solving the momentum-space integral equations [29]. Among scarce studies 
of the low-energy three-body scattering in 2D, the only available calculation in Ref. 25^ 
demonstrated smooth dependence of the 2 + 1 scattering length on the interaction range. In 
the present paper, the precise universal values of the three-body binding energies and the 
scattering length for a particle collision off a bound pai r are calculated. 

It has been known for a long time (see, e. g., j^, 42 1) that the one-dimensional problem 
of N identical particles with the zero-range interactions (in this case, Dirac's 5-function) is 
exactly soluble. On the other hand, the method of the present paper, including the derivation 
of the exact expressions for the coefficients of HREs, is equally applicable to the three-body 
problem in ID. Note that the approach based on solution of HREs was used in Ref. j3| to 
discuss low-energy 2 + 1 scattering in ID. The calculation of the ID three-body problem 
provides a good opportunity to test the approach, to check the numerical procedure, and 
to compare the ID and 2D calculations. For these reasons, the main discussion of the 2D 
three-body problem is complemented by a brief treatment of the corresponding ID problem. 
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The paper is organized in the following manner. The next section contains the informa- 
tion on the low-energy two-body scattering in 2D, introduces the boundary condition model, 
and describes the expansion of the wave function in a set of eigenfunctions on the hyper- 
sphere. Analytical results are collected in Section lllll starting with the eigenvalue equation 
which determines the eigenpotentials for HREs. Furthermore, analytical expressions for the 
coupling terms in HREs are derived. On the ground of these results, the asymptotic form 
of the eigenpotentials and coupling terms is obtained and applied to derive implications on 
the asymptotic behaviour of the solution of HREs. The numerical procedure and the results 
of the numerical calculations are described in Section IIV[ and the last section contains a 
summary and a final conclusion. The results for three identical bosons in ID are briefly 
discussed in Appendix. 



II. METHOD 



The present study is aimed at the description of the low-energy properties of three iden- 
tical 2D bosons with the short-range pair-wise interaction in the limit of the zero interaction 
range. The description turns out to be universal, i. e., essentially independent of the details 
of the two-body interaction. In the low-energy limit under consideration, only the zero total 
angular momentum L = should be considered and only the s-wave two-body interaction 
should be taken into account. The two-body input for the three-body problem is set as 
the universal low-energy description of the two-body interaction by a single parameter, for 
which the two-body scattering length a can be suitably chosen. The scattering length in 2D 
is defined by the asymptotic form of the zero-energy wave function at large inter-particle 
separation r beyond the interaction range, * ~ ln£ This is in analogy with the 

definition of the scattering length in 3D as the distance at which the asymptotic expres- 



sion of the wave function crosses zero, 
the effective-range expansion [l^, 13, 14 
determined by the 2D scattering length a, 

J2i/irk 
fo(k)~ v 



he s-wave scattering amplitude, in accord with 
45^ . in the low-energy limit k — > is completely 



ka ti 1 ' 

In h 7 — i— 

2 1 2 



(1) 



cot So(k) — i 

where k is the wave-number, 5o(k) is the s-wave scattering phase shift, and 7 m 0.5772 is 
the Euler constant. 
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A. Boundary condition model 



In the low-energy limit under consideration, a convenient one-parameter description of 
the two-body interactions is obtained within the framework of the BCM if the interaction 
range is allowed to shrink to zero. The two-body interaction introduced in this way is 
known as the zero-range potential |3| and the Fermi pseudo-potential The equivalent 



approach is also used in the momentum-space representation 



29, 361. Within the framework 



of the BCM corresponding to the vanishing interaction range, the exact scattering amplitude 

fo(k) is determined by the low-energy expression ((J) for an arbitrary k and the two-body 

Ah 2 

binding energy equals -e 27 , which corresponds to the pure imaginary pole of fo(k) at 

ma 1 

ka = i2e~ 1 . Explicitly, the s-wave boundary condition which provides the above-discussed 
low-energy behaviour can be written Q] as 

1 * = . (2) 



lim — 



+o ydr rln(r/a) 

The total interaction of three particles is a sum of two-body potentials, which are replaced 
in the BCM by the two-body boundary condition of the form (J2J) for each pair of particles. 
As the only parameter describing the two-body interactions is the scattering length a, the 
units h — m — a — 1 will be used throughout the paper, thereby the three-body problem 
becomes parameterless. The total wave function \1/ satisfies the boundary conditions and 
the Helmholtz equation, 

[A x + A y + E] * = , (3) 



where x, y is an arbitrary pair of the scaled Jacobi coordinates defined via the particles' 
radius- vectors rj as Xj = rj — and y, = — y= (2r* — — r^). Different sets of the Jacobi 
coordinates are related by Xj = — cxj + syj and y^ = — sxj — cyj, where c = 1/2, s = ±v3/2, 
and the ± sign is chosen if {ijk} is an even or odd permutation of {123}. The wave function 
$ of three identical particles is symmetrical under any permutation of the particles, therefore, 
it is sufficient to impose just one boundary condition, 



lim 

z->0 



d 1 



dx x In x 

where x is any of three inter-particle distances. 



* = , (4) 
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B. Hyper-radial expansion 



Solution of a system of HREs provides an efficien t ap proach to treat both the eigenvalue 
and scattering problem for the three-body system 23|, |35|. This approach is particularly 
advantageous due to the use of the BCM since all the terms of HREs are expressed in 
the analytical form, which allows one to obtain the exact asymptotic form of the wave 
function and to improve the accuracy of the numerical calculations. The system of HREs is 
obtained by expanding the total wave function in a set of eigenfunctions on the hypersphere 

$n(<X,9,R), 

oo 

* = e- R Y,fn(R)®n(a,e,R) , (5) 

n=l 

where the hyper-spherical variables p (0 < p < oo), (0 < < tt/2), and 9 ri (0 < < 7r) 
are introduced by the relations Xi = psina^, iji = pcosaj, and cos^ = (xjVj) / 'xiUi and R = 
In p is a convenient variable in 2D. Different sets of the hyper-spherical variables are related 
by cos 2«j = — c cos 2a j + s sin 2a j cos 9j and sin 2ai cos 9\ = ±s cos 2a j — c sin 2a j cos 9j . By 
definition, $ n (a, 9, R) are regular solutions of the eigenvalue problem on the hypersphere, 
i. e., at fixed R, deduced from Eqs. (J3J) and (jH) 



A 2 + £(i*) 
d 



$ n {a,9,R) = , 



where 



lim 



da a(R + \na) 



$ n {a,9,R) = 



d 2 n n d 
— T + 2cot 2a— + 
a a 1 a a 



2 



sin 2 2a d9 2 



(6) 
(7) 

(8) 



Like the total wave function, the functions $ n (o;, 9, R) are symmetrical under any permu- 
tation of particles, i. e., $ n (a, 9, R) are independent of the index enumerating the Jacobi 
coordinates. 

For each value of the variable R, the problem (jOJ), (J7J) determines an infinite number of 
discrete eigenvalues £ 2 (-R) and corresponding eigenfunctions $ n normalized by the condition 
(<& n |<& m ) = 5 nm . Henceforth the notation (-|-) means integration over the invariant volume 
on the hypersphere dfl = ^ sin 2a da dcos9, where the arbitrarily chosen factor 1/12 is 
suitable for the derivation of the coupling terms in Section IIII Bl The expansion (jSJ) of the 
total wave function leads to a system of HREs which can be written in two equivalent forms, 

d 



lR 2 ~ m) lR 



dR 



Q(R) + V(R) + P(R) - Ee 



2H 



f(R) = 



(9) 



7 



^ + Q(R)) +XJ(R)-Ee 2R 



f(R) = , (10) 



where i{R) is the vector-function composed of the hyper-radial channel functions f n {R) and 
the matrices of eigenpotentials U(R) and coupling terms Q(-R) and P(-R) are defined by their 
matrix elements 

U nm (R) = &(R)6 nm , (11) 

Q nm (R) = ($ n |$'J, (12) 
P nm (R) = (& n \#J » (13) 
with the prime denoting the derivative over R. The identity 

oo 

Pnm ^ ] QnkQmk (1^0 
k=l 

provides the equivalence of the infinite systems of equations in the forms and (jlOJl . 

Although two infinite systems of HREs © and (fTUjl are equivalent, the truncated ones give 
rise to different results, which allows one to estimate convergence with increasing number 
of HREs iV in practical calculations. Notice that N HREs of the form Q reduce to the 
form if the iV-dimensional matrix P^) is replaced by a product of iV-dimensional 
matrices Q^ N \ P^- 1 — > — Q^Q^. It is important that the solution of the truncated 
system of N HREs taken in the form Q gives the upper bound E\ N ^ for the exact energy of 
the zth state E j an d the upper bound for the exact scattering length A, i.e., E[ > Ei 
and AW > A0,Q The proof can be obtained by observing that the truncated system 
of HREs in the form (j2J) can be obtained by application of the variational principle with 
the trial function containing a finite sum of the form (jHJ). On the other hand, the solution 
of HREs (jlOjl . at least in the one-channel approximation, gives the lower bound for the 
ground-state energy j^T | . Solution of the system Q generally provides faster convergence 
with increasing number of equations, while solution of the system (jlOj) does not require 
elaborate calculation of P nm (R). Notice that the scattering length can be calculated by 
solving only the truncated system because the first-channel effective potential U\ (p) 
decreases as 1/p 4 (see Section fill CJ) . In contrast to that, the first-channel effective potential 
in the truncated system (fTTH) is of the form U\ (p) = £f (p) + J2n=i Qin(p) I P 2 an d contains 
a long-range term ~ — l/3p 2 for any finite N, which prevents calculation of the scattering 
length. 
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III. ANALYTICAL RESULTS 



A. Eigenvalue problem on the hypersphere 

It is convenient to take account of the permutation symmetry and to satisfy the boundary 
condition (JJJ) by means of the Faddeev-like decomposition, 

<S>(a,9,R)=j2x(ai,R) , (15) 
i=i 

provided the function x{ a ii R) is symmetrical under the permutation of the particles j and 
k and satisfies the same equation on a hypersphere (JBJ) as the eigenfunction $>(a,9,R). 
The representation (JT5J) is advantageous due to a simple structure of the function x(cn,R), 
which is singular only at one point a = and does not depend on 9 because of the s-wave 
boundary condition. Following Eq. (JJJ), the boundary condition for the function x{°ktR) 
takes the form, 

dx{ai,R) 1 



lim 



doti oa(R + In ocj) ^ ^ 



, (16) 



where the sum contains two functions x{ a j)R) (f° r 3 *)> which are regular in the limit 
ai — > 0. The solution to the eigenvalue problem on the hypersphere is straightforward in 
terms of the Legendre function P v {x) regular at x — 1 

X(a, R) = A(R) P £( R)-i (- cos 2a) , (17) 



where A(R) is the normalization constant. Substituting (JT7j) into the boundary condi- 

the Legendre function as a — > 0, cos 2a) — > 
4^, and calculating the limit cos 2aj^ —> —1/2 as 



tion (Unj), using the asymptotic form o: 
■ sin ttu [In a + 7 + ^(z/ + 1)] + cos 7rz/ 



7T 



«i — >■ 0, one comes to the eigenvalue equation, 

R - 7 - V + I tan |e + ^ec Q = , (18) 

where ip(x) is the digamma function. The same eigenvalue equation, in slightly different 
notation, was derived in Ref. in the limit of the zero interaction range. 

Considering the solution of Eq. (J18|) . it is worthwhile to note that the left-hand side is 



an even function o: 



i.e., R is a function of £ . Similar to the corresponding eigenvalue 



equation in 3D (3J, |35|, the transcendental equation (|18|) determines the infinitely multi- 
valued function £, 2 (R) for an arbitrary complex- valued variable R. In particular, different 



branches of this unique function for the real- valued R form a set of the real- valued £ 2 (-R) 
which play the role of eigenpotentials in the HREs. Hereafter it is convenient to enumerate 
CnC^O by an index n = 1, 2, 3, . . . in ascending order. As R increases from — oo to oo, all 
the terms £ 2 (-R) decrease monotonically in the intervals — oo < £ 2 (-R) < 1, 1 < £f < 25, 
and (2n - l) 2 < £ 2 (i2) < (2n + l) 2 for n > 2. Note that at the exceptional point £ = 3 
the solution of the eigenvalue equation (|18|) gives a finite limit i?o ~ 1-64; nevertheless, 
calculation of the function anc ^ ^ s derivatives in the vicinity of this point requires a 

special care to take into account exact cancellation of divergent terms. 



B. Derivation of the coupling terms 

An important advantage of the BCM is the analytical expression (J18|) for the eigenpoten- 
tials £ 2 (-R) in HREs that allows one to study the asymptotic properties of the solution and 
to simplify the numerical calculation thus improving its accuracy. Evidently, the analytical 
expressions are strongly desirable for the coupling terms Q nm (R) and P nm (R). Whereas 
the direct evaluation of Q nm (R) and P nm (R) by means of the definitions (fT2"j). (fT3*j) is quite 
involved, fortunately, one can circumvent this obstacle by using the explicit dependence on 
the parameter R in the eigenvalue problem (|B|L (|ZJ). Thus, within the framework of the 
BCM one derives the analytical expression for Q(p) and P(p) via eigenpotentials £ 2 (-R) and 
their derivatives over R. 

To simplify the notation, the eigenvalue problem on the hypersphere (JHJ), (J7|) is written 

as 

(A 2 - = , (19) 

lim f^_^ =0 , (20) 
«i^o I dot; a 1 



where e n = — £ 2 + 1, and the function 

M*,R) = ^§ (21) 

ma + R 

tends, as a —>■ 0, to the finite limit which does not depend on the index enumerating the 
different sets of the Jacobi coordinates. Taking the derivatives of Eqs. (fT§j) and (|2*0*|) with 
respect to R, one obtains that $>' n satisfy the inhomogeneous equation 

(A 2 - £ n )$n — £ n®n (22) 
10 



and the boundary condition 



lim ' n 







For derivation of Q n m{R), one starts with the Hellmann-Feynman-type relation 
($ m |A 2 |$;) - «|A 2 |$ m ) = e'J nm + (e m - e n )Q nm , 



(23) 



(24) 



which is obtained by projecting Eqs. (jTUj) and (J22|) onto the functions <&' m and $ m , respec- 
tively. On the other hand, the integrals over the hypersphere on the left-hand side of Eq. (|24|) 
reduce to the contour integrals around three points = in which the functions $ n and 
$^ have singularities. Allowing the length of the contours to shrink to zero and taking into 
account that all three singular points a, = make equal contributions for the symmetry 
reason, one obtains 



($ m \A 2 \K) ~ «|Al*m) = lima 

a— +0 



n da da 



(25) 



Combining the boundary conditions (|2T)J) . (f2~3|) with Eq. (|2*T]) . one finds 



lim a 



$ 11121 _ $' ^ 



(j> n {0,R)<t> m {0,R) 



' n da ~ n da 
and eventually comes from (|2~4*j) - (j2l)j) to the basic relation 

e' n 5 nm + (£m - £n)Qnm ~ </>n(0, R)(f) m (0, R) = . 

The diagonal part of (|77|) provides a simple relation between </> n (0, R) and e' n , 



(26) 



(27) 



(0,i?) = 0, 



(28) 



while the non-diagonal part of ()27)1 combined with ()28)1 gives finally the desired result 

V n m 



(29) 



In a similar way, to derive P nm (R) for n ^ m, one calculates the difference ($^J A 2 |$^) - 
($'J A 2 |$' m ) by projecting Eq. (J22|l onto the functions $' nm and integrating on the hyper- 
sphere, which gives 



(e n - e m )P nm + (e' n + e' m )Q nm = - lim a 



d& d& 



9a 



9a 



(30) 
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In view of fl2D|) and (J23J), the limit on the right-hand side of Eq. flU)) equals m (O, R)(/J n (0, R)- 
0„(O, R)(j)' m (0, R), which allows one to obtain, by expressing <p n (0, R), 4>' n (0, R) via e' n , e" n from 
Eq. (J2BJ), 



(31) 



For derivation of the diagonal terms P nn (R), one requires the functions which satisfy 
the inhomogeneous equation 



(A 2 -£ n )< = 24< + £'>n 



and the boundary condition 



fd$" 
lim IT* 







(32) 



(33) 



Repeating the above procedure to calculate the difference (^'JA 2 ^") — (^"lA 2 !^) and 
taking into account the identity P nn = — one obtains 



V' P 



lim a 

a~*0 



n da 



da 



2[0 / n (O,i?)] 2 -0 n (O, J R)0^(O,/2) , 



(34) 



which after simple algebra in combination with (J28j) gives rise to 

2 



P„ 



QeL 4 I el 



(35) 



The derivation of all the terms in HREs is accomplished with the exact expressions (J29|) . 
(I3T]) . and (j33|) for the coupling terms Q n m{R) and P nm (R) and the eigenvalue equation (JTHJ) 
for £ 2 (-R). Whereas the explicit value of n (O,i?) is of no interest for determination of 
Qnm(R) and P nm (R), it is easy to calculate the limit a — > in Eq. (}2*T]) 

2 



7T 



which, in view of Eq. 



^ n {0, R) = -A n cos ^ n , (36) 

7T Z 

allows the normalization constant to be additionally determined, 
A n = ^ y /-2£ n & sec . (37) 



One should emphasize generality of the derived analytical expressions (|2TJ|) . (JHH), and (|3*K|) . 
The derivation is based essentially on the BCM used to describe the pair-wise interaction. 
Within the framework of the BCM the described procedure is applicable to derivation of the 
coupling terms for a variety of three-body systems in the configuration space of an arbitrary 
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dimension including particles of different masses and scattering lengths and particles with 
internal degrees of freedom. In particular, the analytical expressions of the same form are 
valid for three identical bosons in 3D j2| and in ID (discussed in the Appendix) and for 
three two-species fermions in 3D [soj ] . 

C. Asymptotic expansions and boundary conditions for HREs 

Asymptotic expansions for all the terms of HREs are of interest for qualitative study of 
the described three-body system. In addition, the explicit asymptotic form allows one to for- 
mulate the boundary conditions and to improve the accuracy of the numerical calculations. 
The analytical expressions derived in the preceding sections provide a straightforward deter- 
mination of the eigenpotentials and the coupling terms in the asymptotic region \R\ — ► oo. 

The expansion of eigenpotentials £^(-R) at \R\ — > oo follows from the expansion of the 
eigenvalue equation (jTSjl at the singular points, i. e., near the odd integer £ (except £ = 3) 
and at infinite £. In particular, the expansion at £ — > ioo provides the lowest eigenpotential 
at R — > oo 



The neighboring branches of the multivalued function £(R) are continuously connected at 
infinity so that £ n (R) at R —>■ oo is continuation of £ n _i(i2) at R — > — oo. Thus, the same 
asymptotic expansion at R — > oo for £ n (R) and £ n _i(— R) is obtained by using the expansion 
of .R(£) near the odd integer £, 



As £ n (R) (n > 2) are of the smoothed-step form with the steepest descent at R ~ Inn, the 
asymptotic expansion ()39|) is not uniform in n, viz., it is valid only if R 3> Inn, which hinders 
any consideration of the infinite n limit. Therefore, one needs the asymptotic expansion at 
R — > oo which reproduces the step-like dependence of £ n (i?) at least in the large-n limit 
thus being applicable for both large R and n. The expansion is constructed by using the 
requirement that both the £ n (R) and their derivatives over R for n > 2 coincide with the 



CliR) = -4e 2 ^-^ - - - - e - 2 ^ + 0(e- AR ) . 
via / 3 45 



(38) 




+ 0(\R\- 3 ) , n > 2 . 

(39) 



13 



exact result at the point R n = 7 + ij){n + 1/2) — (— l) n 7rP n ^ 1 / 2 (l/2), viz., one requires 

L 

, which leads 



£ n (Rn) = 2n and g n (R n ) = -4 [vr 2 - 2^'(n + 1/2) + (-l) n 27r 
to the result, 



dv 



u=n-l/2 



TT 



arctana;,, + (— l) n arcsin 



<( J R n )P n _ 1/2 (l/2) 



xl + 1 



(40) 



where x r . 



7T 



-^(i2 n )[i2-7-V(n + l/2)]. As follows from flUJ), £ n (P) (properly 



shifted along both coordinate axes) at large n converge to the function, £ n (R) ~ 2n — 
2 2 

— arctan — (R — Inn — 7). The quite slow (as n -1 / 2 ) large-n convergence is entirely deter- 
mined by the asymptotic behaviour of the Legendre function as v —>■ 00, P J/ _ 1 / 2 (l/2) ~ 



z/ -1 / 2 cos7r(z//3 — 1/4) 49]. Actually, the terms of order ~ n^ 1 ^ 2 contain the dependence 
on n via the expressions (— l) n cos7r(n/3 — 1/4) and (— l) ra sin7r(n/3 — 1/4), which are the 
periodic functions of n with period 3. Thus, one concludes that £ n (i£) up to the leading 
order terms in n belong to three families for different n mod 3. Convergence to the unique 
function is illustrated in Fig. \I\foY two families of £ n (R). 



c 



C 

UJ 1 




-2 2 
R - ln(n) -y 



FIG. 1: Convergence of eigenvalues £, n (R) to the limiting function (bold line). Two families of 
£ n (R) are plotted by solid lines for n = 3m and by dashed lines for n = 3m + 1, m = 1, 5, 9, 15, 25. 



Substituting the above expansions for £ n (R) in the analytical expressions derived in Sec- 



tion UTTBl one obtains asymptotic expansions of the coupling terms. A separate expression 
for the first-channel diagonal coupling term at R — > 00 follows from (J38|) . 

P U (R) = 1/3 + 2/45e~ 4(/? - 7) + 0(e~ 6R ) . (41) 
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Furthermore, using the expansions one finds that Q n m(R) decrease as \R\~ 2 and P nm (R) 
decrease as |-R|~ 4 except the terms Q n i(R) and P n \{R) at positive R, which decrease as 
Q n i(R) ~ P n i(R) ~ e^ R R^ 1 at R — > oo provided R ^> Inn. As discussed above, this 
asymptotic dependence is not uniform in n and one would use the expression (J4(Jj) to obtain 
the uniform expansion which is valid for large n. For example, the desired expansion for 
Qnm(R) at R — > oo takes the form 



4n 



UR) L (-l) n 2x n P n _ 1/2 (l/2) 



1/2 



, ( 2x 
cosh 



-\ -i 



<'n{Rn) 



■ (42) 



Similar to £„(-R), both Q nm (R) and P nm (R) (properly scaled and shifted along the co- 
ordinate axis) converge at large n to the universal limiting functions so that Q n \{R) — > 
(2n)~ 1 / 2 Q l (R - Inn - 7 ), P nl (R) -> {2n)- 1 ' 2 P l {R - Inn - 7 ), P nn (R) -> P(i? - Inn - 7 ), 
and Qnm(R) — > Q(n/m, R — In n — 7) for n > m, where 

Q 1 (y) = (y 2 + n 2 /4y 1/ \coshyy 1 , (43) 

A(y) = Qi(y)(tanhy- 2 ^ 2/ , l , (44) 



Z/ 2 + vr 2 /4 / 



^ 2 



P( y ) = (45) 

V ; 12(y 2 + 7r 2 /4) 2 V ; 

Q(z,y) = z 1 ' 2 {(z 2 - 1) (y 2 + 7r 2 /A) [(y + Inz) 2 + 7r 2 /4]}" 1/2 . (46) 

Splitting of eigenpotentials into three families depending on n mod 3 entails corresponding 
splitting of the coupling terms. The splitting and convergence to the universal limiting 
functions for Q n i(R) and P n i(R) are illustrated in Fig. El 

At infinite separation of particles, i. e., when R — > 00, the leading term of the asymptotic 
expansion for the first eigenpotential is related to the dimensionless energy of the two-body 
bound state so that (, 2 e~ 2R = £ 2 /p 2 — > — 4e~ 27 ~ —1.261, while for other eigenpotentials the 
leading terms in the upper channels are related to the kinematic barriers, ^ 2 e~ 2R —> 1/p 2 
and £, 2 e~ 2R — ► (2n — l) 2 /p 2 for n > 2. Thus, in the asymptotic region, the first-channel 
component of the total wave function describes the two-cluster 2 + 1 configuration, whereas 
the upper-channel components describe the three-cluster configuration. 

Using the expansion of £i(-R) at R — > —00 (JHHJ), one obtains the asymptotic form of the 
first-channel radial function at a small hyper-radius 

MR) ~ e R (R + In £f (r + In ~ - |) . (47) 
15 





R - ln(n) - y 



O 



-2 2 

R - ln(n) - y 




FIG. 2: Two families of Q n i(R) ( a ) and P n \{R) (b) are plotted by solid lines for n = 3m and by 
dashed lines for n = 3m+l, m = 1, 5, 9, 15, 25. The corresponding limiting functions Q\{R— Inn— 7) 
and P\ (R — In n — 7) are plotted by a bold line. In panel (c) the exact result (solid lines) and the 
asymptotic expression (|l2*)) (dashed lines) are compared for the third family of Q n i(R) (n = 3m + 2, 
m = 1,5,9,15,25). 



Given the expansion of £i(R), Eqs. (JTSJ), (|T7|). (J2ZL and the expansion of the Legendre 
function at v — > 0, P v {— cos 2a;) ~ 1 + 2u In sin a 49], the asymptotic form of the first 
eigenfunction on the hypersphere at R — > —00 is 



$i(a,i?) ~ LR + ln 



1 



4 3\ / 4 
i2 + ln- + -J I + In - + ^ In sin «i 



(48) 



As the first-channel contribution dominates in the series (0), the expressions ([47)1 . PH|) entail 
the asymptotic form of the total wave function at R — > — 00, i. e., near the triple-collision 
point, 

' ' x2 / ' X ' Lr i' r ^ r '» ( 49 ) 



4x 2 

* ~ ( R + In - 



f y In sin a, + -R + In - I = In 2 -pin 



3p 2 



In addition, the non-singularity of the lowest eigenpotential in the limit of a small hyper- 
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radius R — > — oo 



exist in 2D 



3, E 



eads to the well-known conclusion that neither Efimov nor Thomas effects 



23). 



For the eigenvalue problem, i. e., for calculation of the bound-state energies, the so- 
lutions satisfy the requirement of the square integrability of the total wave function, 

/oo 
f 2 (R)e 2R dR = 1, and in practice one can use the boundary conditions for the chan- 

nel functions of the form f n (R) — > at R — > ±oo. The asymptotic boundary conditions for 

the low-energy scattering of the third particle off the bound pair are similar to those for the 

two-body scattering. Below the three-body threshold, the wave function in the asymptotic 

region tends to a product of the two-body bound-state wave function <£>(x) and the function 

F(r), which depends on the inter-cluster distance r = V3y/2 and describes relative motion 

of the third particle and a bound pair. At the threshold, i. e., at the zero kinetic energy of 

colliding particles, the 2 + 1 scattering length A is defined by the two-cluster asymptotic 

form F(r) ~ \n(r/A), which leads to the expression \l/(x, y) ~ (p(x)1n — — at y — > oo. 

Taking into account that the first-channel eigenfunction <&i(a,9, R) at a large hyper-radius 

reduces to ip(x)e R and y ~ p = e R for x <C y, one finds the asymptotic form of the channel 

function fi(R), 

2A 

fi{R)~ki-=-R, R^oo. (50) 
v3 

In addition to the asymptotic expression (|5U|) at R — > oo, the first channel function j\(R) — > 

at R —>■ — oo, while all other channel functions f n {R) — ► 0, n > 2 at both limits R — > ±oo. 

The asymptotic form of £f(R) and Pn(R) at R — > oo are of fundamental importance in 

the analysis of the low-energy 2 + 1 scattering. In the lowest channel of HREs, the leading 

term of £f{R) (|38|) cancels the term Ee 2R for the threshold energy E = — 4e~ 27 and the 

next-order constant 1/3 terms of £ 2 (R) (jSHJ) and Pn(R) (jSj) cancel each other, therefore, 

2 

the effective interaction takes the form ^g e_2< ' R_7 ' ) ' wm ch corresponds to the polarization 
interaction V p = — a/(2r 4 ), where r is the distance between a dimer and the third particle 
and a = e 27 /20 ~ 0.1586. This long-range polarization tail of the effective interaction is a 
specific 2D feature (compare, e. g., the exponential fall-off of the lowest effective interaction 
at large distances for three bosons in 3D |35j|). The 2 + 1 scattering length in 2D exists even 
if the effective interaction contains the polarization tail 0, Q (in fact, for the potentials 
decreasing faster than 1/ p 2+s ). This can be seen from the asymptotic solution of the first- 
channel HRE of the system Q at the threshold energy, E = — 4e~ 27 . Up to terms of order 
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0{e 2R ), the first-channel HRE takes the form, 



d2 2 ,-2(fl- 7 ) 



+ — e" 



/i(i2) = 0, (51) 



JR 2 45 

the general solution of (J5"T|) is the linear combination of the Bessel functions 



A(R) ~ d Jo (^2/45 e-^J + C 2 Yq (^2/45 e +7 J • ( 52 ) 

7T C 1 

The asymptotic expansion of the solution (J52j) at i? — > oo, fi(R) ~ — — + 27 — - In 90 — i?, 

2 C2 2 

is of the form (|5U|) . which proves the existence of the scattering length A. As a consequence, 
the leading-order terms of the effective-range expansion for 2 + 1 scattering are of the usual 

7T 

form (Q) for two-body scattering in 2D, viz., — cot 5(k) ~ \n(kA/2) + 7, where k is the wave 
number for the relative motion of a dimer and the third particle. Nevertheless, the higher 
terms of the effective-range expansion are modified by the polarization tail of the effective 
interaction as is known to be the case in 3D scattering [^j]. 

The role of the long-range term ~ e~ 2R or ~ p~ 4 in the first-channel HRE requires 
a careful treatment because there is no clear reason for appearance of the polarization 
potential between a particle and a bound pair. In this respect, it is necessary to study 
a contribution of the upper channels to the effective dimer-particle interactions at long 
distances. Coupling with the upper channels produces in the first channel the nonlocal 

effective potential U C (R, R') which can be estimated in the lowest order of perturbation 

00 d d 

theory as U C {R,R') = £ F n (R)g n (R- R')F n (R'), where F n (R) = Q nl (R) — + —Q nl (R) + 

P n i(R) and g n (R — R') is Green's function in the nth channel. Taking into account that 
(*(R) ~4n 2 , Q nl (R) = (2n)- 1 / 2 Q 1 (R-\nn- 1 ) } and P nX (R) = (2n)- 1 / 2 P 1 (R-\nn -7) for 
large n, one can estimate g n (x) ~ (4n) _1 e~ 2n ' a: ' and F n (R) ~ n~ l l 2 F(R — Inn), where F{x) 

is expressed via Qi(x) and P\{x). As g n (x) — > (2n)~ 2 5(x) for n 00, these estimates entail 

00 

the following local limit of U C (R,R'), viz., U C (R) ~ ^n~ z F 2 {R - In n). Summing over n, 

n 

one finds that the leading term of the effective potential is U C (R) ~ e~ 2R , in other words, 
coupling with the upper channels produces in the first channel the long-range term of the 
same order ~ e~ 2R or ~ p~ 4 as the above-discussed polarization tail. Thus, any conclusion 
on the long-range behaviour of the wave function or, equivalently, on the next-to-leading 
terms of the effective-range expansion for 2 + 1 scattering must be based on the study of a 
large number of HREs. 
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IV. NUMERICAL CALCULATIONS 



The eigenpotentials and the coupling terms P nm (R) and Q n m(R) in HREs were 

calculated by solving the transcendental eigenvalue equation (jTSjl and by using Eqs. ()29|). (J3TJ) 
and ([So]). The derivatives with respect to i? and were replaced by the derivatives 

of the inverse function (dR/d£, d 2 R/d!; 2 , and d 3 R/dl; 3 ) which are easily calculable from the 
eigenvalue equation (|18|). The most involved numerical problem is to calculate the Legendre 
function and its derivatives with respect to the index entering into Eqs. (fTHj). (J2HJ), (JHU 
and (|35j). This is done for both real and imaginary £ by using the Mehler-Dirichlet integral 
representation 



P 



2 



V2 



tt/3 



7T 



(it cos |t 



cost 



(53) 

v _„ 1/2 

for the Legendre function and using for its derivatives the corresponding integral represen- 
tations obtained by differentiating Eq. ()53j) with respect to £. The terms containing an 
integrable square-root singularity are subtracted from the integrand and calculated exactly 
to improve the accuracy. As a result, the Legendre function was calculated with a relative 
accuracy about 1CT 11 whereas the accuracy degraded about one order for each of the sub- 
sequent derivatives. As mentioned in Section 1111 CI accuracy of the numerical calculation 




FIG. 3: The lowest scaled eigenpotentials Q(R)e 2R . Notice different scales for the positive and 
negative R. The arrow marks the two-body bound-state energy — 4e -27 . 

suffers from the subtraction of divergent terms in the vicinity of the exceptional point £ = 3. 
For this reason, ^{R), Q2n(R), and P 2n (R) in a narrow region around the point i? ~ 1-64 
(which corresponds to £2(^0) = 3) were obtained by the interpolation procedure. Under 
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the described approximations, the overall relative accuracy was not worse than 10 -11 for 
the eigenpotentials and 10~ 8 for the coupling terms. It is worthwhile to mention that less 
accurate calculation of the coupling terms is in accordance with a smaller contribution of 
these terms to the final values. The sum rule (|14j) for the coupling terms was numerically 
checked and it was found that the difference Y^k=i QnkQmk — Pnm decreases as N~ 2 with 
increasing N. The eigenpotentials and all the coupling terms for the four lowest channels of 
HREs are shown in Figs. H3 El 
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FIG. 4: Coupling terms Q m n(R) (a), P nn (R) (b), and P nm (R) for n ^ m (c). The arrow marks 
the large hyper-radius limit 1/3 of P±i(R). 



For numerical solution, the truncated system of N HREs is reduced to the form without 
the first derivatives by the transformation f(R) = T(R)f(R), where the orthogonal matrix 
T(i?) satisfies the equation 

dT 



dR 



+ QT = . 



(54) 



Furthermore, one introduces the antisymmetric matrix B by the Cayley transform, B 
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(T — 1)(T + 1) 1 , and solves the equation 

2^| = (B-1)Q(B + 1) . (55) 

This form is preferable because one can use only the upper triangle of the matrices B and Q 
in the numerical calculations, which gives the antisymmetric matrix B and the orthogonal 
matrix T = (1 — B)(l + B)^ 1 independently of the round-off error. Note that in the two- 
channel approximation the non-zero matrix elements of B are explicitly expressed via the 
quadrature, B 2 i = — B\ 2 = tan | / Qi 2 (R)dR. 

Following the described procedure, the truncated system of iV HREs in two forms (JHJ) 
and (fTUj) was numerically solved on the finite interval [R m in, Rmax]- At the first step, 
Eq. ()55|) was integrated and the matrix T = (1 — B) _1 (l + B) was determined at the 
mesh points on [R min , R max \. An arbitrary antisymmetric matrix B serves as the initial 
condition for the matrix equation (J55)) imposed at R m in- The consistency of the numerical 
procedure was additionally shown by checking the stability of the calculated values for 
different choices of the initial matrix T(R min ). Given the calculated transformation matrix 
T(i?), two eigenenergies and the scattering length were calculated by solving the eigenvalue 
problem and the scattering problem at the threshold energy E = — 4e~ 27 for the transformed 
HREs. The zero boundary conditions are imposed in the upper channels, i. e., f n (R m in) = 
f n {Rmax) = for n > 2, whereas the left-end boundary condition in the first channel 
was determined from the asymptotic form of the fi(R) at R — > — oo (@7j). At the right 
boundary, one uses fi(R max ) = for the eigenvalue problem and the asymptotic form (fo^j) 
for the scattering problem. In the latter case, the scattering length is determined via the 

7T C\ 1 

coefficients C\$ calculated at R m ax , viz., In A = — —+27 — In 120, thus taking into account 

2 C2 2 

the polarization tail of the effective interaction beyond the integration region. The boundary 
conditions for the vector- function f(R) were obtained by applying the transformation T(R) 
at the points R m in and R ma x- 

The overall accuracy of the numerical procedure is estimated to provide the calculation 
of the binding energies and the scattering length with the relative error about 3 • 1CT 8 and 
1 • 10~ 6 , respectively. In particular, a sufficient accuracy of numerical integration of HREs 
was obtained by taking R m in = — 14 and R ma x = 1-5, 3.5, and 6.0 for the ground-state, 
the excited-state, and the scattering-length calculations, respectively. The structure of the 
calculated wave function is illustrated in Fig. 03 where the four lowest channel functions f n (R) 
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for the ground state, excited state, and the scattering state are shown. For convenience, the 
solution of the scattering problem is normalized to match the first-channel functions of the 
excited and scattering states at the point R ~ —3.1 corresponding to the first maximum. 
The numerical solution of the truncated system of N HREs provides a set of binding energies 
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FIG. 5: Radial functions of the four lowest channels fi{R) (a), f2(R) (b), fs(R) (c), and fn{R) 
(d) for the ground state (solid lines), the excited state (dashed lines), and the scattering state 
calculated at the two-body threshold energy (dotted lines). For convenience, the radial functions 
for the excited state are multiplied by a factor 5, while those for the scattering state are scaled to 
match at the first maximum the first-channel functions of the excited state and scattering state. 
Linear asymptotic dependence of the first-channel scattering solution is shown by a thin straight 
line in panel (a). 

and scattering lengths, which are presented in Table U in comparison with the calculations 
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40( . It is clearly seen that highly accurate results can be obtained by means of the few- 
channel calculation of the form Q. The contribution to the binding energies from the upper 
channels (for iV > 16) turns out to be comparable with the numerical accuracy. The role 
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TABLE I: The three-body binding energies £o an d E\ (in units of the two-body binding energy) 
and the logarithm of the 2 + 1 scattering length A for identical bosons in ID. The number of HREs 
is denoted by N and the superscripts U and L mark the results obtained by solving HREs of the 
form Q and IjlUj) . respectively. Shown are also the results of fitting the dependence on TV for the 
calculated binding energies and scattering length and those of other calculations. 



N 












; ? 
£\ 








1 A) 

In A 


1 


16, 


.5194096 


16. 


f T c~i o T r\ T 

.5788727 


1 


.26667318 


1 


.29214773 





.891305 


2 


16. 


C~\ t J~\ AAA 

.5219444 


16. 


.5482471 


1 


.26998847 


1 


.27658964 


0. 


f C~\ C~\ i~\ C^i 

.858228 


3 


16. 


.5226064 


16. 


.5302069 


1 


.27033831 


1 


.27263368 





.853238 


4 


16. 


.5226348 


16. 


.5287316 


1 


.27036317 


1 


.27217992 





.851835 


5 


16. 


.5226618 


16. 


.5267981 


1 


.27039042 


1 


.27147416 





.849801 


6 


16. 


.5226787 


16. 


.5249848 


1 


.27040205 


1 


.27101864 


0. 


.848804 


7 


16. 


.5226811 


16. 


.5246296 


1 


.27040405 


1 


.27091797 


0. 


.848343 


8 


16. 


.5226835 


16. 


.5241930 


1 


.27040625 


1 


.27077981 


0, 


.847726 


9 


16. 


.5226854 


16. 


.5237285 


1 


.27040762 


1 


.27066543 


0. 


.847341 


10 


16. 


.5226859 


16. 


,5235979 


1 


.27040796 


1 


.27063093 


0. 


.847125 


12 


16. 


.5226867 


16. 


.5232644 


1 


.27040864 


1 


.27054373 





.846651 


14 


16. 


.5226870 


16. 


5231314 


1 


.27040883 


1 


.27050912 





.846376 


16 


16. 


.5226871 


16. 


,5230155 


1 


.27040895 


1 


.27048209 


0. 


.846186 


oo 


16. 


.5226874 






1 


.27040911 






0.8451 


Ref. [8" 




16.1 ±0.2 




1.25 ±0.05 







Ref. [40] 16.52 1.267 

Ref. [29] 16.522688 1.2704091 

Ref. [25] - - « 1.1 

of the upper channels can be estimated by fitting to the simple power dependence on N, 
which is routinely used in the variational calculations. In the present calculations, it is 
reasonable to fit separately each of three families, i. e., to take into account the periodic 
dependence on N for N = 3m, 3m + 1, 3m + 2. The calculated binding energies are fairly 
well fitted to the a + b/N c dependence for each family with the fitted value of power c ~ 4. 
The logarithm of the scattering length In A converges slower with increasing N than the 
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binding energies, which is manifested by the smaller fitted power c ~ 1 — 1.3. The fitted 
binding energies and the scattering length corresponding to iV = oo are presented in Table H] 
with the overall fitting error in the last digit. As expected, the solution of the truncated 
system (jlUj) provides slower convergence with increasing number of channels N than those 
of the form (JOJ). The calculation based on the solution of the truncated system (jTOjl gives 
a set of binding energies converging as N~ 2 , which is connected with the corresponding 
convergence rate of T,h=i QnkQmk to P nm . 

The calculated binding energies coincide within the declared accuracy with the solution of 
the momentum-space integral equations (3] , which underlines equivalence of quite different 
approaches. The binding energies of a limited accuracy obtained by solving the system 
of HREs Ref. are in agreement with the one-channel calculation of the present paper. 
The older results of js| obtained by solving the integral equations are of low accuracy; in 
addition, the ground-state energy of |8J is above the upper bound found in the present paper. 
The calculations of the 2 + 1 scattering length are rarely available in the literature. The 
present calculation of the 2 + 1 scattering length in the universal limit could be compared 
with the results of Ref. by analyzing the dependence of a% on a 2 shown in Fig. 1 of that 
paper. The three-boson scattering length a 3 is related to the scattering length A defined in 
the present paper as 03 = (2/ir) ln(2A/v / 3), whereas the two-body scattering length a 2 is 
defined in |25j so that the universal limit corresponds to a 2 — ► 0. Considering the smallest 
a>2 ~ 1CT 3 presented as the leftmost point in Fig. 1 of Ref. one obtains S3 ~ 0.8, i.e., 
In A pa 1.1, which is well above the upper bound In A pa 0.8451 calculated in the present 
paper. The discrepancy is presumably because the result of Ref. j3] is not close enough to 
the universal limit and this points to the strong dependence «3(a 2 ) at a 2 — > 0. 



V. SUMMARY AND DISCUSSION 



Universal description of three identical spinless bosons in 2D at low energy is expected 
by analogy with low-energy properties of two particles, which are universal (irrespective of 
a particular shape of the short-range potential) and parameterless if the only significant 
parameter, e. g., the two-body scattering length a, is chosen as a scale. The two-body input 
completely determines the solution near the triple-collision point in the limit of the zero- 
range interactions, therefore, contrary to the corresponding problem in 3D, an additional 
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regularization parameter is not necessary and there are neither Thomas nor Efimov effect 
in 2D. For this reason, a completely universal parameterless description exists in the low- 
energy limit and both three-body binding energies and the 2 + 1 scattering length are the 
universal constants to be determined. 

The BCM is used to describe the pair-wise short-range interaction in the zero-range limit. 
The total wave function is expanded in a set of eigenfunctions on the hypersphere, which 
leads to a system of coupled HREs. The important point is that the analytical expressions 
are derived for all the terms of HREs, which allows one to study the asymptotic behaviour 
and to improve the accuracy of the numerical calculations. One should emphasize that the 
derivation is essentially based on application of the BCM and the Hellmann- Feynman-typ e 
relations, the latter are known to be useful in the calculation of the coupling terms [52]. 
Moreover, the derivation is generally applicable to a variety of three-body problems in ar- 
bitrary dimensions, in particular, the analytical expressions of the same form are obtained 
for three identical bosons j^J and for three two-component fermions in 3D. All the 
considerations and the approach used are equally applicable to description of three identical 
ID bosons, for which the exact solution is known. For these reasons, a brief discussion and 
numerical calculations for the ID case are presented in the Appendix to make comparison 
with 2D results and to check the numerical procedure. 

The analytic expressions are used to analyze all the terms of HREs in the asymptotic 
region, thus obtaining the asymptotic form of the total wave function both for large and 
small inter-particle separation. In this respect, the universal dependence f!49jl is obtained 
for the total wave function in the vicinity of the triple-collision point with the leading 
term ~ ln 3 p and the inter-particle correlations given by In £1X2X3. The large- R asymptotic 
expansions are not uniform in channel number n, therefore, the explicit dependence on n is 
deduced, which reveals convergence of the eigenpotentials and coupling terms to the limiting 
functions of R — Inn = \n(p/n) at large n. The convergence is rather slow and the next-order 
term (~ n~ l l 2 ) in the large-n expansion is periodic in n with period 3; this is displayed by 
observing three families of eigenpotentials and coupling terms, namely, for different n mod 
3. The asymptotic dependence on n is used to study the effect of the channel coupling and 
to shed light on the convergence of the results with increasing number of HREs. One of the 
reasons for slow convergence is the long-range polarization tail ~ e~ 2R ~ p~ 4 of the first- 
channel effective potential and the same order long-range term which arises due to coupling 
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with the upper channels. As a result, one needs to take into account a large number of 
HREs to study the long-range behaviour of the wave function and the next-to-leading terms 
of the low-energy effective-range expansion for 2 + 1 scattering. 

The universal constants, viz., the ground-state and excited-state three-body binding en- 
ergies and the 2 + 1 scattering length, are calculated with high precision by the numerical 
solution of HREs. The binding energies are in excellent agreement (within the declared 
accuracy) with those obtained in the momentum-space calculations j29|, which underlines 
equivalence of two essentially distinct models. The low-energy scattering of the dimer off 
the third particle is completely described by the precise 2 + 1 scattering length. 

In summary, universal low-energy properties of three identical two-dimensional bosons 
are considered within the framework of the BCM used to describe two-body interactions. 
The approach used is based on the solution of a system of HREs, all the terms of which 
are derived in the analytical form. The derivation is quite general and can be applied to 
a number of problems, especially if the interaction is described within the framework of 
the BCM. The asymptotic form of the solutions of HREs is obtained, which allows one to 
describe the wave function both at large and small inter-particle separations. The binding 
energies and the 2 + 1 scattering length of high precision are numerically calculated. 



APPENDIX A: THREE ONE-DIMENSIONAL PARTICLES 



In this appendix, the three-body problem in ID is considered to demonstrate general 
applicability of the approach used, to check the numerical accuracy, and to compare conver- 
gence of the ID and 2D calculations. The choice is based on well-known exact solubility of 
the one-dimensional N-body problem with the zero-range interactions j^l, 42 1). As usual, the 
problem becomes parameterless by introducing the natural units h = m = 1 and by choosing 
the potential strength to fix at unit values both the two-body binding energy, e% = 1, and 
the two-body scattering length, a — 1. The exact result for the binding energy of n identical 
particles in ID 4l|, |42| is e n = |n(n 2 — 1). One should also mention that the ground-state 
wave function of three identical particles is of a simple form ty gs = C exp(— Y^k \ x k\), where 
the scaled Jacobi coordinates and are introduced similar to the above-discussed 2D 
case. The solution at the threshold energy E — — 1 determines the wave function of three 
particles \l/ sc = J2k exp(— \xk\) — 4exp(— | \ x k\), which entails infiniteness of the 2 + 1 
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scattering length or existence of the zero-energy virtual state 

Thereafter, the approach described in the paper is applied to calculate the three-body 
binding energy €3 and the 2 + 1 scattering length A of three identical particles in ID. The 
wave function satisfies either the equation 



d 2 d 2 3 



dx 2 



i=l 



m = 



(Al) 



where the zero-range interaction is a sum of the Dirac 5-functions, or the free equation 
complemented by the boundary condition that can be written for the each pair of the 
identical particles as 



lim 

x^±0 



dx 



^ = . (A2) 
Similar to Section Hi Bl one introduces the variables p and aj and expands the wave function 



P 



-1/2 



£ f n (p)$ n (a, p) 



(A3) 



n=l 



in a set of eigenfunctions $ n (a,p) on a circle of constant p, which leads to the systems of 
ordinary differential equations for the functions f n {p) which are analogous to Eqs. Q, (fTUj). 
Eigenpotentials in these systems are defined by the solution of the eigenvalue problem on a 
circle and the coupling terms are defined by the analytical expressions of the same form (|29|) , 
(J31|) . and (J35j) as in 2D, provided the derivatives are taken over p. Recall that the derivation 
of the analytical expressions for the coupling terms in Section fill Bl is equally applicable in 
ID. 

For the symmetry reasons, the eigenvalue problem on a circle can be solved in the interval 

< < vr/6 by imposing the zero boundary condition — — = 0ata = 7r/6 and the 

da 

boundary condition at a = 0, 



lim 



d 

da ^ 



^ = 



(A4) 



which follows from Eq. (|A2|) . The solutions of the eigenvalue problem on a circle satisfying 
( d 2 \ 

the equation I ^— ^ +£^J $ n (a,p) = take a simple form <3> n (a,p) = B n cos(a — 7r/6)£ n , 
where the eigenvalues £ n (p) are defined by the transcendental equation, 

. (A5) 



7T 



£ + pcot -f 
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Due to simple dependence £(p) (|A5|) . one can derive simple analytical expressions for the 
coupling terms, for example, 



cos 4 § 
p = - 

1 nn 



x (x 2 — 3) (x — 2 sin x) — 6x 2 cos x + 3 sin 2 x 



3x 2 (x + sinx) 4 



(A6) 



7f 

where x = — £ n (p), as follows from Eq. ([3*5|). 

Similar to the 2D case, the numerical solution of the HREs in the form Q, (fTUjl with 
zero boundary conditions gives the three-body binding energy, whereas the solution at the 
threshold energy E = — 1 gives the 2 + 1 scattering length. In the latter case the asymptotic 
form of the wave function is a product of the two-body bound-state wave function (fi2{x) = 
exp(— |x|) and the function F = 1 — -^4, which determines the asymptotic form of the 
first-channel function fi(p) = p _1 / 2 (l — 2 a) - ^ s snown i n Table ITT1 the calculated binding 
energy rapidly converges to the exact value £3 = 4 with increasing N, whereas the calculated 
scattering length rapidly grows with N, which manifests infiniteness of the exact scattering 
length. Both and are fairly well fitted to the a + b/N c dependence with the fitted values 

TABLE II: The three-body binding energy €3 (in units of the two-body binding energy) and the 
2 + 1 scattering length A for identical bosons in ID. A number of HREs is denoted by N and 
the superscripts U and L mark the results obtained by solving HREs of the form Q and (jlOj) . 
respectively. 

N e 3 ^ e 3 L A 

1 3.99934308 4.00728928 3.32000 • 10 2 

2 3.99998993 4.00055763 7.9633 • 10 3 

3 3.99999902 4.00013463 4.555 • 10 4 
5 3.99999994 4.00002429 3.31 • 10 5 
7 3.99999999 4.00000816 1.1 • 10 6 
9 4.00000000 4.00000367 2.6 • 10 6 
12 4.00000000 4.00000149 7.9 • 10 6 
15 4.00000000 4.00000074 > 10 7 

of power c ~ 6 and c ~ 4, respectively. The fitting of the scattering-length dependence on 
shows that the calculated A grows as A^ 3 . A better precision of the ID calculation in 
comparison with the 2D one is basically due to a simple form (jA4|) of the eigenvalue equation 
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that provides a better accuracy of the eigenpotentials. Both for ID and 2D calculations, 
the second type of truncation of the HREs provides energies converging to the exact values 
from below. 
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